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Abstract. We discuss an asymmetric energy model (AEM) introduced by Giardina et al. in [10] . 
This model is expected to belong to the KPZ class. We obtain lower bounds for the diffusion 
coefficient. In particular, the diffusion coefficient is diverging in dimension one and two as it is 
expected in the KPZ picture. 

1. Introduction 

In their well-known paper ([H]), M. Kardar, G. Parisi and Y.-C. Zhang introduce a model for 
the evolution of the profile of a growing interface 

d t h = V 2 h+-(Vh) 2 + r J (x,t) (1) 

Here h(x,t) S M is the height of the interface at location x G M and at time t > and rj(x,t) 
is a space-time white noise. Starting from a flat state at time t = 0, they are interested in the 
evolution of the fluctuations 

f(s,t) =< [h(x,t)~ < h(x,t) >] 2 > 1/2 

The intuitive picture is that the width grows with time as a power law up to a saturation time that 
scales with the substrate size L as L z where z is the scaling exponent. In other words, we expect 
that 

t(L,t) = L 2 -*f(t/L*) 

where the scaling function f(x) saturates at large x and f(x) ~ x( 2_2; )/ z for x ~ ([10]). By 
dynamic renormalization- group techniques, Kardar et al. show that in dimension 1, the dynamic 
scaling exponent is z = 3/2 and z = 2 in dimension d > 3. The dimension d = 2 is the critical one 
and their numerical studies indicate 1/z = 0.620.04. As noticed by Kardar et al., equation (pQ) can 
be mapped to the Burgers equation for a vorticity-free velocity field 

d t v + v ■ Vv = V 2 v - Vt](x, t) (2) 

with v = — V/i and ij(x, t) a space-time white noise. 

Burgers equation is also closely connected to driven diffusive systems. Consider a diffusive 
system under constant uniform driving force described by a nonlinear Langevin equation. Then 
a quadratic order expansion of the density gives the Burgers equation Q. In [3], H. van Beijern 
et al. investigate the steady-state scattering function for driven diffusive systems with a single 
conserved density. Mode-coupling arguments predict that in dimension d = 1 (resp. d = 2), density 
fluctuations spread as t 2 / 3 (resp. t(logt) 2 ' 3 ) whereas they are of order t 1 / 2 (ordinary diffusion law) 
in dimension d > 3. Guided by ideas of universality, we expect that a large class of microscopic 
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models whose evolutions are in a suitable coarse time and length scale an approximation of the KPZ 
equation dU) (or the noisy Burgers equation) have a universal scaling exponent z and a universal 
scaling limit. 

Asymmetric Simple Exclusion Process (ASEP) is a natural discretization (see |20j ) of the sto- 
chastic Burgers equation. The dynamics are given by asymmetric random walks on Z d with a drift 
in some direction such that jumps of particles to occupied sites are forbidden (exclusion rule). It 
can be reinterpreted as a growth model which is a natural discretization of the KPZ equation. 
During the last decade, a lot of work has been accomplished to test the validity of universality 
predictions for ASEP. In the one dimensional case the value of the dynamical exponent z = 3/2 has 
been confirmed. Not only the exponent but also the scaling function was obtained ([?])■ Moreover 
the limit is the same as the largest eigenvalue distribution in the random matrix theory (see [5] and 
references therein for more informations). Nevertheless an important fact has to be mentioned: the 
methods used to obtain dynamical scaling exponent and limit distribution are very dependent of the 
specific properties of ASEP. Indeed the main results are valid and proved in the one dimensional 
case and for the Totally Asymmetric Simple Exclusion Process (TASEP), which corresponds to 
nearest neighbors jumps in the right direction. These results can not be carried for general ASEP. 
For example M. Prahofer and H. Spohn (|_18j) compute the current fluctuations for the TASEP but 
their proof does not work for other ASEP models. A more robust method has been introduced 
by M. Balazs and T. Seppalainen ([2]) but is restricted to attractive systems. Even if they do not 
obtain the scaling limit function, they are able to establish the order of current for the nearest 
neighbors ASEP but not for general ASEP (obtained in [20] by generalized duality techniques). 
Recently this method has been developed in the context of the Asymmetric Zero Range Process 
(AZRP) (see [I]). 

To my knowledge, the only models belonging to the KPZ class and for which a rigorous proof 
of the scaling order has been obtained are ASEP, AZRP, Poly Nuclear Growth model (PNG) and 
related models ([2], [I], [6], [TO]). In dimension 2, the class is even more restrictive and the only 
rigorous result is the scaling order obtained by H.T. Yau for general ASEP f[23|). Hence, the class 
of microscopic models for which one can rigorously prove they belong to the KPZ class is very small 
and it is hence of extreme importance to have simple models for which one can rigorously prove 
they are in the KPZ class. 

The aim of this paper is to study a non- attractive model introduced by Giardina et al. in 
|10| and to show it presents anomalous behavior in low dimension as it is expected in the KPZ 
picture. In [10] . Giardina et al. consider Symmetric Energy Model (heat conduction model in their 
terminology) and show the system has a dual process. 

They also introduce an asymmetric generalization of the model that we call the Asymmetric 
Energy Model (AEM). AEM should belong to the KPZ universality class. It presents several 
analogies with ASEP but also differences ([TO]). In this paper, we develop generalized duality 
properties for AEM and obtain lower bounds for the bulk diffusion coefficient D(t), i.e. the variance 
of the two points correlation function (see ©). KPZ approach predicts large time behavior of D(t). 
In particular, D{t) is expected to be divergent in dimension 1 and 2 and finite in dimension d > 3. 
The goal of this article is to obtain lower bounds for D(t) consistent with this (theorem 1). The 
proof of this result is based on generalized duality techniques introduced by C. Landim and H.T. 
Yau in the context of ASEP. They have been developed in several directions but essentially for 
lattice gas dynamics. Our main sources of inspiration are given by [I] and [15j . 

The paper is organized as follows. In section 2, we define AEM and introduce the diffusion 
coefficient D(t). Section 3 is devoted to the generalized duality properties of the process. Section 4 
contains the technical lemmas necessary for the proof of theorem 1. The proofs of the main theorem 
are given in section 5,6,7 for the 1, 2 and d > 3 case. The paper is ended by remarks in section 8. 
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2. The Asymmetric Energy Model (AEM) 

The system is composed of atoms indexed by x S The canonical basis of M. d is denoted by 
(ei, . . . , ed)- Each atom has a momentum p x £ R. Momenta are exchanged during the stochastic 
evolution in such a way that kinetic energy is conserved. The generator L = S + A of AEM is 
defined by 

d 

»=1 

and 

d 

(Af)(p) = J^Ci ^ PsJfc+ei -p x d Px+ei ) 

Here p = (px) x ez d i s an element of the state space = M z and / is a smooth local function of p. 
Parameters (aj)i=i ... d regulate the strength of the asymmetry in each direction. Let 

Mdp) = ® xeZd (2vrT)- 1 / 2 exp(-pl/2T)dp x 

be the Gaussian product measure with temperature T. [it is an invariant probability measure for 
L. Moreover S is symmetric and A is antisymmetric in L 2 (f2,/iT). We fix now T > and denote 
fir by < • >. Energy of site x is denoted by E x = p 2 . The formal total energy ^2 x& %d E x is a 
conserved quantity of the dynamics and one has 

d 

L(E X ) = 2A(E X ) + 2Y^a l V £t (E x _ et E x ) 

i=l 

where A is the usual (f-dimensional discrete Laplacian: 

d 

AE X = {Ex+ei + E x - ei — 2E X } 
i=i 

and V ei is the discrete gradient in the direction a: 

\/ g - ^Ej jp £ / g ^ — £ — |— g ■ £ / £ / £ / j° g ■ 

The microscopic instantaneous current in the direction ej is given by 

jx,x+ ei {p) = 2Ve.(E x ) + 2(iiE x E x+ei 
and one has the following microscopic continuity equation 

d ft 

E x (t)-E x (0) = y ^ / (V e 

i3x—ej,x 

) (p(s))ds + M x {t) 

i=i ^ 

where M x (t) is a martingale. 

We are interested in the energy-energy correlation function S(x,t) defined by 

S(x,t) =<E x (t)E o (0) > -T 2 

At equilibrium, the mean value of the current in the direction e« is ji(T) = fiT{jx,x+ei) = 2a{F 2 . 
By the conservation law (cf. [22], pp. 263-264), one has the two following informations for the 
average location and velocity of the structure function: 



5^S(x,t) = 2T 2 



and 

— 2^xS*(x,t) = tv 



where v = Ya=i fi( T ) e i = ^ t Ya=i a i e i- 

The third natural quantity to study is the bulk diffusion coefficient which is defined by 
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= ^rAYl ^S(x,t) - 2T 2 (vit)( Vj t) \ (3) 
Based on mode coupling theory (|_3J), it is expected that 

d=l 

D(t) ~ I (logt) 2 / 3 , d = 2 
1, d > 3 

for large i. Let Wi(p) be the normalized current in the direction 

w<(p) = Jo,ei(p) - k{T) - j'i{T){pl - T) 

or more explicitly 

Wi(p) = 2(1 + ai T)V et (E ) + 2 ai (E - T){E ei - T) 
For local functions / and g in L 2 (f2), we define the semi-inner product 

< /,0 >= (< ^/<? >-</><<?>) 

Here is the usual shift on $7. The semi-norm corresponding to <C ■, • 3> is denoted by || • ||. Note 
that (discrete) gradient terms g = T x h — h and constants vanish in this norm. 

A formal integration by parts gives the following formula for the diffusion coefficient (see [13]) 



DiAt) = %+ ^ « f m(p(s))ds , /%(p(s))ds » (4) 



2 4T 2 

The Laplace transform of the diffusion coefficient is then given by 

poo -J. i 

^ e- A *tAj(t)^ = ^2 + ^^2<C^,(A-L)- 1 ^» (5) 

We have to mention that we don't have a rigorous proof of equality (j3J. In fact even the existence 
of the equilibrium infinite volume dynamics is a non trivial problem ([8]). Several methods exist in 
the literature ( [9] , [T7] and references therein) but they are not directly applicable and we plan to 
extend these methods for AEM in a future work. In the rest of the paper, we assume the following: 

Assumption 1. The operator L defined on the set of local integrable smooth functions of 0, is 
closable and its closure also denoted by L is the generator of a strong Markov process. Moreover 
the set of local smooth integrable functions on Q is a core for L. 

In the following RHS of (gj) will be used as definition of D(t). For ASEP, the validity of (gj) 
can be established by coupling techniques one can note translate in our context (see |14j for more 
informations on this subject). It follows that in terms of Laplace transform, behavior of D{t) for 
large t is, in a Tauberian sense, equivalent to behavior for small A of 

< Wi, (A - LY 1 w i > 



For simplicity, we will restrict ourselves to the case i = j. We will prove the following theorem 
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Theorem 1. There exists a constant C > such that 

'c\- x i\ d=l 

<C Wi, (A - L) _1 u>; >> < C\ logA| 1/2 , d = 2 

C, d>3 

Moreover, if d > 3, we /iaue also an upper bound 

< u>;, (A - >< C' 1 

In a Tauberian sense, this theorem means that D(t) > Ct 1 / 4 for d = 1, D M (t) > (logt) 1 / 2 for 
d = 2, C" 1 > > C > for d > 3. 

Observe that the assumptions are not so relevant when we are only interested in lower and upper 
bounds. Without any assumption one can obtain similar bounds if we define the diffusion coefficient 
by a finite volume limit procedure. It means we can define D(t) by the following limit (when it 
exists): 



XiXjS(x,t) - 2T 2 (vit)( Vj t) 



where is the d-dimensional discrete torus of length N and the dynamics is now defined on ~R t n . 
Then one can prove the following lower bounds 

f00 , , fCA-V* d=l 

A 2 / e- xt tDrp(t)dt > < C\ logA) 1 ^ d = 2 
[C, d>3 

where C > is independent of A and N. Moreover we have the corresponding upper bound if 
d > 3. The advantage to deal directly with the infinite volume definition of the diffusion coefficient 
is that it simplifies notations and avoid to work with discrete Fourier transform but directly with 
continuous Fourier transform. 

3. Duality 

For simplicity, we fix the temperature T equal to one and we denote by fx or < • > the standard 
Gaussian product measure The Hilbert space associated to is denoted by L 2 (S1) and the 
corresponding inner product is denoted by < -,- >. Let h n ,n > be the sequence of Hilbert 
polynomials which is an orthogonal basis in L 2 (R, (2ir)~ 1 l 2 e~ x / 2 dx): 

p -x 2 /2 



J 



h n (x)h m (x) 



= 5 mn n\, n,m>0 



The dual space of AEM is N and elements £ of N are seen as configurations of the " generalized 
dual process" . If £ = (£ x ; x G 1< d ) is a configuration of the dual space, we will say that £ is local if 
£ x / only for a finite number of sites x £ Z d . In such a case the number of particles |£| is defined 
by 

Kl = E & 

x& d 

For any local configuration £ : Z d — > N, we define the multivariate Hermite polynomial function 
by 



Hdp) = n 



where n : N — > N\{0} is a suitable normalization function we will precise later. Remark that 



r<5. 



Here, the notations £! and n(£) are for IXcfeO anc ^ rL n (£*0- 

Any local function / G L 2 (f2) can be decomposed uniquely as a finite linear combination of local 
functions 

£ n>0|£|=n 

The coefficients of this linear combination are given by a real valued function F defined on the 
dual space N z . Such a function is said to be of degree n > if F(g) = as soon as |£| / n. For 
example, the normalized current Wi(p) has the following decomposition 

Wj = 2aj <5g. + gradient term (6) 

where & is the configuration with two particles on site and two particles on site ef. 

'2 if x = 
= S 2 if x = e, 
otherwise 

V 

We know examine how acts the generator on the dual space. For this purpose we need to introduce 
notations. We introduce a cemetery configuration o which does not belong to the dual space. Any 
function F : N z —> K is extended to a function F on N z U {o} by F{6) = 0. If £ is a configuration 
belonging to the dual space N^ d and x,y are two sites of Z d then £ x >+ 2 >f> -2 is the configuration 
obtained from £ by moving two particles from site y to site x. If £ y < 2 then ^ x < +2 <y~ 2 = o. The 
configuration £_ x > +2 (resp. i x '~ 2 ) is the configuration obtained from £ by adding (resp. by removing) 
two particles on site x. In the second case, if £ x < 2 then £_ x ~ 2 = o. If / £ L 2 (J7) is a smooth local 
function such that 

/ = E = E E F ^ 

then 

d d 

w = J2 + a ^f} = E { E(^ F )(o^ + o< E(^ F )(o^ 
i=i i=i [ ^ a 

We take 

n(jfe) = fc!! = fc(fe — 2)(fc — 4) ... , n(0) = l 
A long but elementary computation (see appendix) shows that 

tern) = Y,& + 1 K*-ed F (e' +2 ' x - ei '- 2 )-F(o] 



+ E^ + 1 )^[^ 



3?,-(-2,X-|-6j 



2 ) - ^(0] 



and 
with 



X, 2,33 6j,-|-2 



)] 



and 

x 

and 



x 

Remark that S = Yj%=i * s the generator of a Markov process reversible w.r.t. the measure 



n 2 (H) 

This Markov process is not irreducible but if we restrict this process to the invariant subspace 

£e,n = U; ^1 = n, & = e x mod 2} , e G {0, 

then the restriction is irreducible. Let TL £ ,n be the subspace of functions F vanishing outside £ e>n . 

Remark that the process corresponding to the generator S is the same as the dual process derived 
in [10]. Nevertheless it is important to observe that our basis {Hg} is different from the basis of 
|10j . The key advantage of our choice is that the basis {H^} is orthogonal and the computations 
simplify considerably. 

The operators A®, Af and A~ are not Markov generators. The operator A® conserves the degree 
and is antisymmetric in L 2 (m). The operator Af increases the degree by 2 and A~ decreases the 
degree by 2. In L 2 (m), one has (A^)* = —Af . 

We will use the following abusive but very convenient notations in the sequel. If /, g £ L 2 (f2) 
are smooth local functions with coefficients in the basis given by local functions F and G then 
we have: 

<f,g>=<F,G> 

with < F,G > defined by 



We write similarly 
with 



<f,g>= n$G(€Mt) 

< /, 9 >=< F, G > 



« F, G »= E F (0G(r x 0m(0 (7) 

x&L d geN zd \{0} 

where t x £ is the shifted configuration by x, meaning (t x £) z = £ x +z- Remark that in ((TJ), the sum 
is carried over configurations £ with at least one particle. 

4. Free particles Approximation 

We introduce the generator A (discrete Laplacian) of independent random walks on 7L d . It is 
given by: 

(AG){rj)= 7] x [G(rf'~ 1,y ' +1 ) - G(r])] 



\x-y\=l 

Here r\ € N z and rj x '~ 1,y,+1 is the configuration obtained from r\ by moving a particle from site x 
to site y. G is a real valued local function defined on N zd . A configuration rj can be seen as an 
element of /S n where S n is the symmetric group of order n. The identification is given by the 
map 

(xi,...,x n ) G (Z d ) n /^ n ^ V 
7 



with 

n 



i=i 

A function G in the domain of A is then identified with an element of ©^ =1 L 2 ((Z d ) n /S n ) (with 
respect to the counting measure). We denote the standard inner product on ®^ =1 L 2 ((Z d ) n /£„) by 
< ■) ' >o: 

VF,GGL 2 ((^)«/S„), <F,G> = Ex=(x 1 ,...,„) 6 (^ W s/WG(x) 

1 (8) 
- ^Ex=(xi,... A )e(z")» f ( x ) G ( x ) 

We define also the scalar product with translations associated to < •, • >o 

« F, G » = (< F > r * G >o " < F >o< G > ) 

xei d 

We recall that A is a positive self-adjoint operator w.r.t. < •, • >o- 
Lemma 1. Let F(£) be a local function belonging to Tto^n and define G{rj) = F(2n) then 

C- lc T 2n < -AG, G >o<< -SF,F>< Cn < -AG, G > 
where C is a positive constant independent of n and F. 
Proof. One has 

<-SF,F>= Y, E U^ + Q[m x '- 2 ' y ' +2 )-F(0] 2 m(0 

?6£o,n \x-y\=l 

For any integer k, we have 
and 

( + } " ((2fc + I)!!) 2 " n} =1 (2j + 1) " l\ V 2, + 1 
so that 2" fe < m(21fe) < 1 and (2/3) fc < m(2k + 1) < 1. It follows that 

2~ 2n < inf m(£) < sup m(£) < 1 

We use now the fact that 1 < (1 + £ y ) < (1 + 2n) and we conclude. □ 
It follows that the same lemma is also true for the inner product with translations since 



Lemma 2. Let F(£) a local function belonging to TCo^n and define G(rj) = F(2rf). There exists a 
constant C independent of n and F such that 

c -i 2 -2n ^ _ AG) G >0 < < _ S F, F »< Cn <C -AG, G > 



The iJ_i x norms are defined by 

\\Ff_ j A =< F, (A - S^F >= sup{2 < F, G » - < G, (A - S)G >} 

G 

and 

llFHij A o =< F, (A - A)~ X F > = sup {2 < F, G > - < G, (A - A)G > } 

G 

In these formulas, function F is a local function from N zd into R and the supremum is carried over 
local functions. 

It follows easily from lemma 2 and equation (|8|) that there exists a positive constant C(n) such 
that for every local square integrable function F belonging to TCo^n 

^yll^ll-i,A,o < ll*1l-i,A < C{n)\\mFf_ 1X0 
In this inequality, the function mF is defined by (mF)(£) = m(£)F(£). 

Every function belonging to to Wo,2n can be seen as a symmetric function F from (Z ) n into R: 

F(xi,...,x n )=F(24 1 + ... + 2^J 

where 5 X denotes the configuration with only one particle on site x. Sum of two configurations and 
multiplication by an integer are defined in the standard way. In the sequel, we will identify F with 
F. 

To obtain a lower bound, we use the following variational formula for the Laplace transform (see 

my- 

< Wi, (A - L)~ l Wi > 

= su P/ {2 « /, Wi » - « /, (A - S)f » - « Af, (A - S)"M/ »} (10) 
= sup F {2 « F,% » -\\F\\l x - PF|| 2 _ 1>A } 

and we restrict the supremum over functions / = F(^)H^ such that F belongs to Wo,4- 

Lemma 3. Let d > 1 and F, G : (Z d ) n — ► R 5e symmetric local functions. There exists a constant 
C := C(d,n) such that 

\< G,A°F >| < G(G, (A-SG) 1/2 (F, (A - S)F) 1/2 
The same inequality is valid for the inner product <C •, • 

|< G,^°F >| < C < G, (A - 5)G > 1/2 « F, (A - S)F > 1/2 
Proof. By definition of m, for any nonnegative integers k, £, we have 

m(k + 2)(k + 2)m(£ - 2)(£ - 2) = (Jfe + l)m{k)lm{£) (11) 
and by simple changes of variables, one has: 

(A1F, G) = + 1) [F(^ x ~ 2 ' x+ei ' +2 ) - F(0] - G(e'- 2 ' x+e - +2 )] m(0 

- £*fe- e * + 1) [F(e'- 2 ^- e - +2 ) - f(o] [g(o - Gr'- 2 ^- 6 ^ 2 )] m(o 

The estimate is then a simple consequence of Schwarz's inequality. The second part of the lemma 
follows from the definition of <C ■, • 3>. □ 

Corollary 1. Let d > 1 and F : {'L d ) n — ► R 5e a symmetric local function. There exists a constant 
C = C(d,n) > suc/i that 

P°F|| 2 < G||F|| 2 



Recall we want to obtain a lower bound for the formula (jlOp . Corrolary Q] implies we can forget 
the term ||«4°i ?, ||^_ 1 A in the variational formula (|10|) . 



5. The one dimensional case 

We recal that every function F belonging to Tio^n is identified with a symmetric function F 
from Z n into R: 

F( Xl ,...,x n ) =F(2S X1 + ... + 26 Xn ) 

In the sequel, we forget the and consider a function F £ Hon as a symmetric function on Z n . 
To obtain a lower bound, we restrict the supremum over functions / = such that F 

belongs to Tto : 4. We identify hence F with a symmetric function on Z 2 . 
By corollary [H ||^4°i ? || 2 _ 1 A is of no importance and can be forgotten. 

Since F £ Wo, 4, we have that A + F belongs to Wo, 6- A + F is identified with a symmetric function 
on Z 3 . 

For all symmetric function F(x,y) (resp. G(x,y, z)) from Z 2 (resp. Z 3 ) to Z 2 (resp. Z 3 ), we 
define 

F*(a) = ^2F(z,z + a) 

resp. 

z, z) 

zez 

One can check that 



«F,F» =^(F») 2 , «G,G» = £ (G*(a,[3)) 2 



and 



« f, (-af) » = E t F > + e ) - F »] 2 ' « G ' (- AG ) >>o= E E t G > + e ) - 

where the sum over e is carried over all e e Z 2 sucht that e = ±(0, 1), ±(1,0), ±(1, 1). 

The advantage of this notation is to reduce the degree of functions. We take now a function 
F G Wo, 4 (i-e. a symmetric function on Z 2 ) and we compute G* = (mA + F)*. One has 

3 

(mA + F)(x, x,x + l) = -^{F(x, x + 1) - F(x, x)) 

3 

(mA + F)(x, x,x-l) = -(F(x, x) - F(x - 1, a?)) 

(mA + F)(x - 1, x, x + 1) = ^(F(x, x + 1) - F(x - 1, a;)) 

(mA + F)(x, x + 1, y) = ^(^(^ + 2/) ~ F ( x > 2/))' y^ar-l,a;,aT±l,a;±2 
(m^4 + F)(x, y, z) = otherwise 
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(13) 



A long but elementary computation shows that: 

G*(x,y) = ^[5 l (y-x)-5 1 (x-y)](F*(y)-F*(x)) 

- ±8 1 (\x\)6 1 (y-x)(F*(x + l)-F*(x)) 

- ^i(|y|)5i(x-y)(F*(y + l)-F*(y)) 
+ l(Si(x)h(v) ~ S^(x)S^(y))(F*(0) - F*(l)) 
+ ^-i(x)6>2(y)(F*(y-l)-F*(y)) 
+ i5^(y)5> 2 (x)(F*(x-l)-F*(x)) 

+ ^x(x)5> 3 (y)(F*(y)-F*(y + l)) 

+ ^i(y)S> 3 (x)(F*(x)-F*(x + l)) 

We choose now the following test function F: 

F*(x) = A-V4 e -A 3 /*|*| (14) 

and evaluate for this test function the terms appearing in the variational formula (|1U|) . In the 
sequel, C denotes a positive constant independent of A which can change from line to line. 
By using the decomposition (JSJ) of Wi, we have that for small A > 0: 

<<5 Sl ,F>= 2ai^F(x,x + l) = 2aiF*(l) ~ A" 1 / 4 

X 

Moreover, since F belongs to Ho,4, there exists a positive constant C such that 

C- 1 < F, F >o<< F,F ><C < F,F > 

and by Q, the same is true for the inner-products <C ■, • 3>o / •C •, • S>. The norm of F with respect 
to <C -, • is easy to evaluate and is of order A -1 / 4 . The third term to estimate is <C F, —SF 
By lemma O this term is of the same order as <C F, —SF 3>o Thanks to formula (|12|) and explicit 
form of F* , we obtain that <C F, (-S)F 3> is of order A 1 / 4 . Hence we proved 

2<%,F>~ A~ 1/4 , A<F,F>~ A~ 1/4 , < F, -SF >~ A 1/4 (15) 

We now evaluate the last term 

< K.4+F)], (A - A) -1 [m(yl + F)] > 

Lemma 4. Let F*(x) = A _1 / 4 e~ A3/4 l a: L We have 

< [m(.A + F)], (A - A)- 4 [m(^ + F)] » ~ A^ 4 (16) 
Proof. For any function / : Z n — ► R, we introduce the Fourier transform / of / defined by 

m = E f ee[o,if 

By (fT2|) . we have 

|G*(«,«)I 2 



« w^i. (a - A)-' [m (^F)] »„= y Ma + 4e(i) + 49(s + () 
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where 6(u) = sin 2 (7™). We have 

« K^F)], <* " A)" »„< ^t J G ^f im = (G-, (A - A)-' G ') 

We can express Fourier transform of G* in terms of F* . We write G* = G\ + G% + ■ ■ ■ + Gg with 
G*, . . . , Gg the eight terms appearing in f)13[) . We claim that 

f« Gj, (A - A) _1 Gj » ~ j = 1,5,6,7,8, 

\« Crj-, (A - A) _1 Gj » = O(AlogA), j = 2,3,4 

We begin by the proof of the first claim. We have 



and 



G{(x,y) = \[S 1 (y - X )-S 1 (x- y)\ (F*(y) - F*(x)) 



&l(s,t) = -^{sin(2vrs) + sm(2irt)} F* (s + t) 



The Fourier transform of F* is easy to compute and is given by 

A-V4(i _ 2e - 2A3/4 ) 



F*(s) 



(l-e- A3/4 ) 2 + 4e- A3/4 sin 2 (7rs) 
and we have 

\F*(s)\ < ^ 

\ P W - A3/2 +S in 2 (7T S ) 

as soon as A is sufficiently small. Moreover, a simple computation shows that 

f i » \Gl(s,t)\ 2 >in(27r S )+sin(27rt)) 2 |F*(s + t)| 2 

hoir dsdt x + 4e(s) + 4e(t) ~ ho,i? dsdt \ + 46( s ) + 4e(t) (if) 

= J^ds\F*(s)\ 2 sin 2 (tts)U x (s) 



cos 2 (7™ I 



= / n — n vr du 

A + 2 sin 2 — + 2 sin 2 ' 



where 



2 / V 2 

We are interested in the behavior of (|17p as A goes to zero. The critical points of the integrand are 
and 1 and by symmetry arguments, we can restric ourselves to 0. By equation (4.12) of [1], one 
has 

C 

U x (s) < 



It follows that 

f 1 ' 2 ~ „ 9 f 1 ' 2 sin 2 (7rs) 

4 / ds\F*(s)\ 2 sm 2 (7rs)U x (s) < CX / ds . V 

Jo Jo (A 3 / 2 + sm (7rs)) 2 v A + s 2 

Standard analysis shows that this last term is of order A -1 / 4 which proves the first claim. 
Let us examine the term Gs (terms G5, Gq, G7 are evaluated in the same way). We have 



G* s (x,y) = A 1/4(1 4 6 A ] 5 l (y)5 x > z e- x3/Ax 
12 



By a direct computation, we obtain 

\G* 8 (s,t)\ 2 <C X 

\ 6 i z + sin [tvs) 

It follows that 

«O 8 ,(A-A)-. Os » = kipMt —^ML— } 

r i 1 f „i dt 

ds 

< cx r! — ; 

= 0(A -1 / 4 ) 

because we have 

nl dt C 

< 



I 

Jo 



, x + d( s ) + e(t) ~ ^x + W) 

)r the terms G2, 63,64, we use 

that 



For the terms 62,63,64, we use the fact that \G*(s,t)\ < 6A 1 / 2 if j = 2,3,4 and we observe 



dsdt 

O(logA) 



W \ + 6(s) + 6(t) 

The lemma is proved. □ 

We can now conclude the proof of theorem [1] for the one dimensional case. We use the variational 
formula (fTUl) and the test function 



Here a is a positive constant we will fix later and the function F is defined in (|14p . By (|15p . lemma 
H] and equation (fTUl) , we obtain 

<. wi, (A - L)-^! >> Cia\-^ 4 - C 2 a 2 X 1/4 - 6 3 a 2 A~ 1/4 

where 61,62,63 are positive constants independent of a and A. If a is chosen sufficiently small, 
the lower bound is of order A -1 / 4 and theorem [1] is proved. 

6. The 2-dimensional case 

For simplicity, we assume 02 = and a\ = 1. 

Lemma 5. Fix R > and assume d > 2. There exists a positive constant 6 = 6(d) independent 
of A and n snc/i £/iai /or any local symmetric function F : (Z d ) n — > E 

n 

E (l|x l -x J |+|x fc -x,|<^ 2 ) o < 6n 3 (F, (A - A)F) 
i,j,k,l=l 

Proof. This lemma is proved in lemma 4.2 of [15J. □ 

Lemma 6. Let d > 2 and n > 1. There exists a constant 6(n) such that for any local symmetric 
functions F : (Z d ) n=Fl -> M and 6 : (Z d ) n -» R 

KC^F) - (0,^)0 1 < C(n) (F, (A - A)F)J /2 (G, (A - A)G)J /2 
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and it follows that 

|< G,A ± F > - < G,A ± F > | < C{n) «F,(A- A)F > 1 / 2 < G, (A - A)G >q /2 
Proof. One has 

n 

k,i=l 

with x fc the vector x where the coordinate x fc has been removed. Observe it is the operator defined 
in formula (3.4) of |15j . We have 

n 

(G,A+F) - (G,AfF) = £ l {xfc=Xf+ei} G(x)(m(x) - l)(F(x fc ) - F(x*)) (18) 

k,e=i 

If x G (Z rf ) n is such that Xj 7^ Xj for every i 7^ j then m(x) — 1 = 0. Hence we can introduce the 
sum of the following indicator function 

in the sum (|18p . We recall that |m(x) — 1| < C(n) for a constant C(n). It follows that 

|<G,.A t +F>- <G,.A+F > | 

< c ip) E El^Wi^-x^iiix^x,^!!^)-^)! 

One concludes by Schwarz inequality. □ 
Corollary 2. Lei d > 2 and F : (Z rf ) ra — > IR. There exists a constant C(n) such that 

P^H^A < C(n) {W^f-XM + H^lllAo) 
Proof. We have only to prove the lemma with the inner product <C • , • 3> replaced with the standard 
inner product without translations < •, • > thanks to ([9]). Recall that 

(A + F, (A - 

= sup {2 < A + F, G> -\<G,G> - <G, (-S)G >} 

G 

Then we use lemma [6j lemma Q] and equation ([8]) to conclude. □ 



We can now complete the proof. Recall that we want obtain a lower bound for the RHS of (|10p 
and remark that in the dual basis, one has 

««>!,/»= 2oi E F ( 2 ^ + 2 ^i) 

where F is the function corresponding to / in the dual basis, we restrict the supremum over 
functions belonging to Wo, 4- Hence, by corrolarydj corrolary [2] and lemma [21 one has: 

<C wx, (A - L) _1 wi » 
> sup J401 Y j F(2z,2z + 2e l )-C\\F\\\ xo ~ C\\A + F\\\ X A 

The maximizer for this new variational problem can be computed using Fourier transform (cf. 
lemma 3.3 of |15|). It turns out that 



<C wi, (A - L) 1 wi >> 



C\\og\\ 1 ' 2 d = 2 
C d>3 
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7. The ^-dimensional case for d > 3 



In this section, we show that the diffusion coefficient is finite and strictly positive if the dimension 
d is greater than 3. For the asymmetric simple exclusion process, this has been proved in [16]. The 
lower bound follows from the previous section. The upper bound is obtained by ignoring the 
asymmetric part of the generator: 

<C wi, (A - L)~ l wi » (19) 

= sup {2 « /, Wl » - « /, (A - S)f » - « Af, (A - S^Af »} 
/ 

< sup {2 « /, wi » - « /, (A - 5)/ »} 
/ 

We write this last variational formula in the dual basis and we recall that 

« Wl J »= 2oi J2 F ( 26 z + 2t 

z&L d 

with f(p) = ^£ F(^)H^{p). Let us decompose F in the orthogonel sum composed of subspaces 
7in :£ (s are elements of {0, l} z<i ): 

F — ^ F n ^ £ , F n<£ £ 7~(-n,e 

Since S send an element of TC n ,e ° n an element of W n ,e ; we have 

ll-^lll.A = H-^n-.elll.A 

It follows that in the third line of (|19|) . one can restrict the supremum over functions F belonging 
to 7^2 0- If F belongs to 7^2,0 then F is identified with a symmetric function on (Z d ) 2 and we have 

ii^Hm ^ c \\ F \\i,xv 

for a positive constant C. The supremum is then easily computed using Fourier transform. In the 
supremum appears the Green function of the discrete Laplacian which is finite only for dimension 
d > 3. We have 

« wi, (A - >< C, d>3 

8. Conclusions and final remarks 

For d = 2, the diffusion coefficient is expected to be of order (logi) 2 / 3 . It has been proved by 
Yau in [23] in a very technical paper for ASEP. It is not clear that the method of [23] can be applied 
for AEM. Indeed, constants C(n) appearing in lemma Q] are exponential in n and the method of 
|23j seems to be restricted to polynomial dependence in n. Remark that ASEP is the only model 
belonging to KPZ class for which such behavior is proved for d = 2. For the case d > 3, it would 
be interesting using generalized duality techniques to establish a fluctuation-dissipation equation 
(|16j) for AEM, meaning a decomposition of the current in the form 

Wi = V^i + Lhi (20) 

for functions (fi and hi in a suitable Hilbert space. 

In this paper we obtained lower (and upper for d > 3) bounds for the diffusion coefficient for 
AEM. The strategy was based on generalized duality techniques similar to [3] and [15]. A recent 
paper of Balazs and Seppalainen (cf. [2]) improves considerably the lower bounds obtained in [15] 
in the one dimensional case for the diffusion coefficient. For nearest neighbors ASEP (but not for 
general ASEP), the authors of [2] are able to prove upper and lower bounds with the right order. 
Nevertheless, their method has restrictions : a key ingredient is attr activity of the process (AEM is 
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not) and the method has only been developed for the one dimensional case. Of course, generalized 
duality techniques have also restrictions. The key ingredients are: 

• If (n p )p, is a shift invariant family of stationary measures indexed by the conserved quantity 
p (e.g. the density), then L 2 (// p ) can be decomposed in an orthogonal sum ®„>oW„. 

• The symmetric part of the generator S maps 7i n into 7i n (S conserves the degree). It means 
that the symmetric part of the generator restricted to TC n is the generator of a reversible 
Markov process with a finite number of particles. 

• In general, the asymmetric part A of the generator does not conserve the degree but the ac- 
tion of A on each subspace H n is a bounded operator from (H n , ||-||i) into (®^=^ no 'Hj, II • 
where no is a fixed positive integer. 

A model introduced by Seppalainen in [2T] belongs to the KPZ universality class and should 
have anomalous behavior in low dimension as ASEP and AEM. Even if one consider this attractive 
process, methods of [2] are difficult to apply (there is not concept of second class particle for this 
process). This process is also difficult to study with duality techniques because the symmetric 
part of the generator does not conserve the degree. But a suitable modification of the process 
can be studied with duality techniques. The slight modification is the discrete counterpart of 
the asymmetric energy model. The symmetric part is given by the KMP process (|12j) and the 
asymmetric part by the asymmetric part of the process defined in |21j . The state space is M zd . 
For a real valued local function /(£) defined on the state space of the process, the action of the 
generator L on / is given by 

(l/xo = x)Er*tf(^^ (21) 

i=l x€% d 

where the exchange operator E and the transfer operator T are defined by 



E x ,y,p£ = £ + \p(£ x + Cy) - 6r]4 + [(1 " p)(U + Q ~ QS y 
T x ,y,p£ = C+P£,x(5y ~ S x ) 

This Markov process is well defined, conserves the energy £ x and centered exponential product 
measures of the form 

X 

are invariant. The symmetric part of the generator is the KMP process. At the difference of 
the symmetric part of the process defined in [2T], it conserves the degree and generalized duality 
technique presented here (the basis is composed of multivariate Laguerre polynomials) can be 
applied to the process defined by (I2ip and one obtains similar lower bounds for the corresponding 
diffusion coefficient. 
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9. Appendix 

We derive here the expression of the generator L in the dual basis H^. Recall that the Hermite 
polynomials (h n ) n satisfy the following equations: 

d?h n dh n 

—r^r - u— h nh n = 

au z au 

dh n _ 

*7 " nh ^ (22) 
K+i(u) = uh n (u) - nh n -i(u) 
h = 1 

where we adopt the convention that h n = if n < 0. We denote by the multivariate Hermite 
polynomial without the normalization factor n(£): 

Hz = n(0Hz 

Fix x ^ y in 7L d . By (f2"2"j) . we have 

(Pydpx ~Pxdp y )H^ = ^ x H^x,-l,y,+ l — ^ y H^y,-l,x, + l 

It follows that 

(Pyd Px - Pxd Py ) 2 Hi- = £ x (£ x - l)H^,-2, y ,+2 - £ x (£ y + 1)H^ 

and we get finally 

(Pyd Px ~ Pxd Py fH^ = - + 2)%,- 2lHl+2 - + 1)^ 

+ (& + 2)(e w - 1)%,+2,„,- 2 - e„(£ x + 1)^ 

because we have — 2)/n(k) = 1/k. To obtain the expression for (<Sj.F)(£), we write / = 
F(£)H£ and we have 

= E F (0{(^ - !)(£»+«< + 2)H^- 2 , x+H ,+2 - Uix +ei + } 

+ {& + - i)%^.^,- 2 - ex+ ei fe + 1)^} 

Recall here that our conventions are such that = if there exists z E Z rf such that £ 2 < 0. By 
a suitable obvious change of variables, we get 

+ E |e^ + [F(e- +2 ' s+e ^ 2 ) - no] | H t 

The computations are similar for the asymmetric part. We have 

PxPy{Pyd Px -Px9p y )H(: = P x P y (£ x H^,-l,y,+ l - £ y H£y,-l, x ,+l] 

By symmetry we can restrict ourselves to compute 

PxPy{ixH^,-i, y ,+i) 
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which is equal to 

Cxflfv.+a + + l)(Cx - l)%,-2 + - l)%,-2, H ,+2 + + 

and we obtain that 

PxPy(p y d Px -p x d Py )H^ 

is equal to 

£ x Il£y,+2 + ix{iy + l)(£r - l)-ff^,-2 
-iyH^x,+2 - iy{ix + l){£y ~ l)H^y,-2 

+& - Zy)Hz 

+£,x(Z,X — l)H^x-2,y,+2 — £y(£y ~ 1)-H^»,-2,b,+2 

Since H% = n(^)H^, we obtain 

PxPy{Pyd Px -p x d Py )H{: 

= ix{i y + 2)Hgu,+2 - £ y (t; x + 2)H^,+2 

+ {ty + !)(£» " l)H^-2 - (Zy - + 1)^,-2 

+ (£,X — + 2)H^x,-2,y, + 2 — (£y — l)(Cx + 2)H^x, + 2,y-2 

If / is a local smooth function such that /(p) = F(^)H^ then 

+ {(£* - l)(6s+e, + 2)%,-2,, +ei ,+2 - (£ x+e . - + 2)%,+2,, +ei ,- 2 } 

By suitable changes of variables, we get the announced expression for A4F. 
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